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^the  usual  pressure  gradients  found  in  engineering  applications.  The  method  includes  a newly 
formed  wake  modification  function  for  the  rimilarity  laws  as  well  as  added  effects  due  to 
low  Reynolds  numbers. 

The  variation  of  momentum  loss  boundary-layer  thickness,  local  skin  friction,  and  local 
velocity  profile  may  be  calculated  for  the  boundary  layers  on  hydrofoils  and  two- 
dimensional  appendages.  The  various  formulations  are  shown  to  correlate  well  with 
available  experimental  data. 
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ABSTRACT 


Analytical  relations  have  been  derived  for  calculating  a developing  two- 
dimensional  turbulent  boundary  layer  in  a pressure  gradient  from  two  simul- 
taneous differential  equations:  momentum  and  shape  parameter.  An 
entrainment  method  is  used  to  obtain  the  shape  parameter  equation.  Both 
equations  incorporate  the  velocity  similarity  laws  that  provide  a two-parameter 
velocity  profile  general  enough  to  include  any  range  of  Reynolds  numbers. 

The  entrainment  factor  is  based  on  the  characteristics  of  equilibrium  pressure 
gradients  in  a way  which  can  accommodate  the  usual  pressure  gradients  found 
in  engineering  applications.  The  method  includes  a newly  formed  wake- 
modification  function  for  the  similarity  laws  as  well  as  added  effects  due  to 
low  Reynolds  numbers. 

The  variation  of  momentum  loss  boundary-layer  thickness,  local  skin 
friction,  and  local  velocity  profile  may  be  calculated  for  the  boundary  layers 
on  hydrofoils  and  two-dimensional  appendages.  The  various  formulations  are 
shown  to  correlate  well  with  available  experimental  data. 

ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  authorized  and  funded  by  the  Independent 
Research  Program  of  the  David  W.  Tayloj  Naval  Ship  Research  and  Development  Center 
under  Project  ZR-023^01 , Work  Unit  1-1 541-002. 

INTRODUCTION 

Analytical  prediction  of  turbulent  boundary  layers  in  pressure  gradients  for  two- 
dimensional  flows  such  as  exist  on  foils  or  struts  may  be  expeditiously  accomplished  by  the 
so-called  integral  methods  that  involve  ordinary  instead  of  partial  differential  equations.  Here 
the  classical  von  Kirmin  momentum  equation  is  usually  accompanied  by  a differential 
equation  for  the  streamwise  variation  of  the  shape  parameter  of  the  velocity  profile.  Such  a 
shape-parameter  differential  equation  may  be  obtained  readily  from  the  entrainment  method 
originally  introduced  by  Head.1 

As  developed  by  Head,  a one-parameter  relation  is  empirically  correlated  between  the 
usual  shape  parameter  and  a newly  defined  entrainment  shape  parameter.  Then  another 
newly  defined  entrainment  factor  is  empirically  correlated  with  the  shape  parameter  in  a 


1 Head,  M.R.,  “Entrainment  in  the  Turbulent  Boundary  Layer,"  Aeronautical  Research  Council  RAM 
(Sep  1958). 

A complete  listing  of  references  are  given  on  page  40. 
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one-parameter  relationship.  A two-parameter  method  was  later  developed  by  Head  and  Patel2 
based  on  the  Thompson  charts  for  a two-parameter  velocity  profile.  The  entrainment  factor 
based  on  equilibrium  pressure  gradients  is  presented  on  graphs. 

A two-parameter  entrainment  method  based  on  velocity  similarity  laws  has  been 
presented  by  Granville3  as  part  of  a synoptic  group  of  such  integral  methods.  The  entrain- 
ment factor  is  obtained  from  boundary  layers  in  equilibrium  pressure  gradients.  By  being 
stated  in  terms  of  momentum  thickness  Reynolds  number  and  shape  parameter,  the  entrain- 
ment factor  is  decoupled  directly  from  the  pressure  gradient  and.  hence,  becomes  more 
responsive  to  nonequilibrium  pressure  gradients. 

Recently  an  improvement4  has  been  developed  for  the  law  of  the  wake  and  thus  in  the 
analytical  description  of  the  similarity  laws.  It  is  now  proposed  to  apply  this  modification  to 
the  entrainment  method.  In  addition,  the  skin  friction  law  to  be  used  is  to  be  based  on  the 
similarity  laws  directly  instead  of  a previously  used  generalized  version  of  the  Ludwieg- 
Tillmann  formula.  This  allows  inclusion  of  roughness  effects  and  drag-reducing  polymer 
effects  more  readily.  When  the  turbulent  boundary  layer  is  considered  a distinct  layer,  a 
momentum  equation  and  an  entrainment  equation  may  be  derived  from  the  governing 
Reynolds  equations  for  turbulent  flow  without  the  usual  thin  boundary-layer  considerations. 
The  results  are  additional  terms  arising  from  cross-pressure  variations  and  normal  stress, 
turbulent  fluctuations.  These,  however,  are  important  usually  only  close  to  separation  of  the 
boundary  layer. 

Velocity  similarity  laws  incorporating  a modification  to  the  law  of  the  wake  are  used  to 
obtain  a iwo-parameter  variation  of  entrainment  shape  parameter  as  a function  of  the  usual 
shape  parameter  and  the  local  Reynolds  number  based  on  momentum  thickness.  The  velocity 
similarity  laws  are  also  used  to  obtain  a variation  of  local  skin  friction  with  these  same  two 
parameters,  even  at  low  values  of  the  Reynolds  number. 

As  stated  before,  the  entrainment  factor  is  analytically  derived  from  boundary  layers  in 
equilibrium  pressure  gradients  from  a two-parameter  relation  with  shape  parameter  and  local 
Reynolds  numbers.  This  decouples  the  entrainment  factors  directly  from  the  local  pressure 
gradient  and  provides  a response  to  the  cumulative  history  of  the  developing  boundary  layer. 


JHead,  M.R.  and  V.C.  Patel,  "Improved  Entrainment  Method  for  Calculating  Turbulent  Boundary  Layer 
Development,”  Aeronautical  Research  Council  RAM  3643  (Mar  1968). 

1 Granville,  P.S.,  "Integral  Methods  for  Turbulent  Boundary  Layers  in  Pressure  Gradients,”  Journal  of  Ship 
Research,  Vol.  16,  No.  3,  pp.  191  -204  (Sep  1972). 

4Granville,  P.S.,  “A  Modified  Law  of  the  Wake  for  Turbulent  Shear  Flow,”  NSRDC  Report  4639  (Sep  1975). 


Hence,  the  entrainment  factor  based  on  equilibrium  pressure  gradients  may  be  applied  to 
nonequilibrium  pressure  gradients. 

For  convenience  all  the  analytical  relations  required  for  calculating  the  development  of 
the  turbulent  boundary  and  the  associated  velocity  profiles  are  listed  in  the  appendix. 

TWO-DIMENSIONAL  TURBULENT  BOUNDARY  LAYERS 
COORDINATE  SYSTEM 

For  the  two-dimensional  boundary  layer  on  foils,  the  coordinate  system  is  defined  in 
terms  of  families  of  surfaces  parallel  to  the  foil  surface  separated  by  normal  distance  y.  The 
other  coordinate  s is  specified  along  the  contour  representing  families  of  surfaces  orthogonal 
to  the  foil  parallel  surfaces.  The  arc  length  s is  measured  in  the  stream  direction  from  the 
forward  stagnation  point  while  the  normal  distance  y is  measured  outward  from  the  foil 
surface. 


EQUATIONS  OF  MOTION 


For  turbulent  flow,  the  Reynolds  equations,  which  are  the  averaged  Navier-Stokes 
equations  for  laminar  flow,  are  applied  to  the  boundary-layer  region.  For  the  foil  coordinate 
system  and  with  the  usual  neglect  of  terms  involving  longitudinal  curvature,  the  Reynolds 
equations  of  motion  are  in  the  s-direction.5 


3u  + 3u  _ 1 8p  | I 3r 

3s  3y  p 3s  p 3y  p 3s 


and  in  the  y-direction 


3v  3v  1 dp  1 dr  I do 

u— +u  — = -—  r1-  + - — - - — 

3s  3y  p 3y  p 3s  p 3y 


The  equation  of  continuity  is 


(2) 


(3) 


5Schlichting,  H.,  “Boundary-Layer  Theory,”  Sixth  Edition,  McGraw-Hill  Book  Co.,  New  York  (1968). 


where  u = velocity  in  s-direction 

v = velocity  in  y-direction 
r = shearing  stress 

o,  = turbulent  normal  stress  in  s-direction 
a = turbulent  normal  stress  in  y-direction 
p = pressure 
p = density  of  fluid 


MOMENTUM  EQUATION 


The  momentum  equation  results  from  an  integration  in  the  y-direction  of  Equation  (I), 
of  the  s-equation  of  motion,  and  the  inclusion  of  Equation  (3),  the  equation  of  continuity 


where 


dfl  „ 9 dU  rw  Id  r6  , . 1 d Cb 

ds  IH  2>  U ds  -^2  jJ0  <»  - 


dy 


rl 

0 = momentum  thickness  = J 

“'o 

5 = boundary-layer  thickness 
H = shape  parameter  = 6 */6 


dy 


6*  = displacement  thickness 


ss‘/S('-u) 


dy 


tw  = wall  shearing  stress 
pg  = pressure  at  outer  edge  of  boundary  layer 
U = streamwise  velocity  at  edge  of  boundary  layer 


(4) 

(5) 

(6) 
(7) 


This  momentum  equation  differs  from  the  classical  von  Karman  equation  by  the  inclusion 
of  the  cross-pressure  and  normal-stress  terms. 

The  momentum  equation  may  also  be  conveniently  written  as 


d 0 
ds 


+ (H  + 2)fj 


(8) 


where 


the  cross-pressure  term 


(9) 


and 


the  normal  stress  term  a * — / ot  dy  . (10) 

ds70 

It  should  be  noted  that  the  cross-pressure  and  normal-stress  terms  become  important  usually 
only  close  to  separation. 

ENTRAINMENT  EQUATION 

The  entrainment  equation  has  been  developed  by  Head1  on  the  basis  of  physical 
arguments  regarding  the  growth  or  entrainment  of  fluid  in  the  developing  turbulent  boundary 
layer.  It  may  also  be  obtained  analytically  by  integrating  the  equation  of  continuity;  see 
Equation  (3). 


/•5|,|9u  sv 

i *dy  0 


which  by  the  Leibnitz  rule,  and  the  boundary  condition  v = 0 at  y * 0,  gives 

d r8  . ,,  dS 

-_/  udy-U-t’s-O 

where  vg  = v at  y = 5,  the  entrainment  velocity. 

To  determine  the  entrainment  velocity  vg,  Equation  (1),  the  equation  of  motion,  is 
evaluated  at  y = 5;  then 


Since  in  genera!  along  a streamline 


du  l 3u  t 8u  dy 
ds  3s  dy  ds 


then  at  y = 6 


Likewise, 


/ 3u\  _ dU  /du\  d6 
\9s /g  ds  \9y/6  ds 


VELOCITY  SIMILARITY  LAWS 

The  velocity  profile  for  turbulent  boundary  layers  in  pressure  gradients  may  be  specified 
by  the  classical  velocity  similarity  laws  which  imply  a two-parameter  system.  Correlation 
with  experimental  data  has  been  well-established  s 

For  the  usual  moderate  pressure  gradient,  the  inner  law  or  law  of  the  wall  of  smooth 
fiat  plates  (zero  pressure  gradient)  has  proved  applicable  which  is  stated  in  nondimensional 
ratios  as 


fly*] 


(23) 


(S' 


where  u = streamwise  velocity  component 

uT  = shear  velocity 

UT=  Jrjp 

tw  = wall  shearing  stress 

, = ^y 


(24) 


(25) 


v = kinematic  viscosity  of  fluid. 

The  inner  law  may  be  extended  to  include  roughness  effects6  and  drag-reducing  effects 
of  polymer  solutions.7 

For  pressure  gradients,  the  outer  law  or  velocity-defect  law  may  be  stated  as 


¥->M 


(26) 


where  B2  is  a velocity-defect  factor  that  varies  mostly  with  pressure  gradient.  Roughness 
effects  or  drag-reducing  polymer  effects  do  not  affect  the  outer  law. 

Within  the  boundary  layer,  the  two  similarity  laws  may  be  considered  to  overlap  which 
by  the  classical  Millikan  analysis  leads  to  logarithmic  functions  within  the  overlapping  region3; 
see  Figure  I. 


— = A In  y*  + B 
ur 
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(27) 


or 


u_  6 * 


7 -Hi), 


(28) 


where  A = constant 

Bj  = the  law-of-the-wall  factor,  which  is  constant  for  smooth  surfaces  and  a 
variable  for  rough  surfaces  and  drag-reducing  polymer  solutions6,7 


^Granville,  P.S.,  “The  Frictional  Resistance  and  Turbulent  Boundary  Layer  of  Rough  Surfaces."  Journal  of 
Ship  Research,  Vol.  2,  No.  3,  pp.  52-74  (Dec  1958). 

7GranviUe,  P.S.,  “The  Frictional  Resistance  and  Velocity  Similarity  Laws  of  Drag-Reducing  Polymer  Solutions,' 
Journal  of  Ship  Research,  Vol.  12,  No.  3,  pp.  201  -212  (Sep  1968). 
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y0*  * inner  limit  of  log  law 
= outer  limit  of  log  law 


The  region  in  the  boundary  layer  outside  the  logarithmic  region  has  been  fitted  by  the 
Coles  law  of  the  wake,  which  in  subsequent  modification4  may  be  included  with  the  log  law 
for  the  range  y*<y*<v 


yo  y 

or  for  the  range  — < -7  < I 
V « 


U - u . , y _ w fyl  fy 

— = -A.n|+  B2-B2I  I -Aql| 


where  w is  a wake  function,  and  q is  a wake-modification  function. 
Moses4  gives  a polynomial  expression  for  w 


and  Granville4  for  q 


-(f)2  H) 


For  convenience  w and  q are  considered  to"  start  at  the  wall  y = 0. 


BOUNDARY -LAYER  PARAMETERS 


Boundary-layer  parameters  suggested  by  the  outer  similarity  law  are 


which  from  definitions  results  in 


"■•(f) 


where  o 


U /'wV,/2  • 
UT  \pU2/ 


is  a local  skip-friction  coefficient. 


which  from  definitions  results  in 


■-TfirHi] 


'.-•’ft'-!) 


The  Rotta-Clauser  shape  parameter  G defmed  by 


becomes  from  previous  definitions 


G = l2/l, 


* ° ” ff) 


This  relation  for  G is  independent  of  the  choice  of  velocity  profile. 

If  the  similarity  laws  as  expressed  by  Equations  (27)  and  (30)  arc  used  to  evaluate  I, 
and  l2.  then 

''  = nAtTB2*7  1 


1,  , i£I2  A2 , 2U  . B . 13  2 , :>i°  ~ *: 

2 2520  140  AB2  35  B2  n 


where 


X,  =(Aln  yo*  + B,  -A)y0*-  Pl  ± 

'/o  ur 


X2  = ((Alny<;  + B1)2-2A(Alny0*  + B1)  + 2A2|y0*-J/"yo^2  dy*  (43) 

For  very  high  values  of  Reynolds  numbers,  tj  becomes  increasingly  larger:  hence,  the 
X-tcrms  may  be  neglected.  Then 


(44) 
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and 


4819  2 
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(45) 


The  same  results  are  obtained  if  the  outer  law  is  assumed  to  hold  all  the  way  to  the  wall 

y - 0. 


ENTRAINMENT  SHAPE  PARAMETER 


The  entrainment  thickness  is  defined  as  the  difference  5-5*  and  the  entrainment  shape 
parameter  H is  defined  as  the  ratio 


For  one-parameter  velocity  profiles, 

H = f[H] 


and  for  two-parameter  velocity  profiles. 

H = flH,  Rfl) 


where 


(46) 


(47) 


ONE-PARAMETER  RELATIONS 


The  power-law  velocity  profile 


gives 


and 


H 


2H 
H - I 


dH  _ 2 

dH  (H-l)2 


(48) 


(49) 


(50) 


II 


An  empirical  fit  to  airfoil  data  for  1.4  x I03  < < 7.6  x 104  by  Head8  yields 

another  one-parameter  relation 


and 


H = 1.535  (H  -0.7)-2-715  + 3.3 


~pj=-4.l68(H-0.7r3-715 


(51) 


(52) 


TWO-PARAMETER  RELATIONS 

The  similarity  laws  provide  a two-parameter  relation  for  the  entrainment  shape  parameter. 
From  Equations  (34)  and  (39) 


and  with  Equation  (39) 


Then  from  Equation  (46) 


BB) 

£ - li  /H  - 

6 G { H / 

MrfHirn)-" 


(53) 


(54) 


(55) 


This  equation  has  been  also  obtained  by  Michel.  Quemard.  and  Durant.® 

G/lj  may  be  related  to  G through  Equations  (44)  and  (45)  with  B,  as  the  implicit 
parameter.  At  separation  00  and  G/lj  -*  1.4857.  A close  fit  is  given  numerically  by 


g2.75 


(56) 


Then  the  entrainment  shape  parameter  becomes 


B-(irn)  ( 


1.4857  + 


0.4739  A 5 A2 


&)- 


(57) 


where  G is  a function  of  H and  through  Equation  (39).  Then. 


* ‘‘Proceedings- Computation  of  Turbulent  Boundary  Layers,  1968  AFOSR-IFP  Stanford  Conference," 
Edited  by  Rline,  SJ.  et  al.,  Vol.  I,  Thermosciencet  Division,  Stanford  University  (1968). 


• * 


and 


3H  _ H(H  - 2)  - H / H2 
3H“  H(H-I)  "VH-I 


/0.4739A  ( 13.75  A2\ 

T 1 .a 

L 3 In  o] 

V G G2-75  / 

H(H-  1) 

Y 3H~  J 

3H  = _ / H2  \ /0.4739  A 13.75  A2\  3 In  a 
In  Rg  \H  — 1/  \ G q2.75  / 3 In 


(59) 


LOCAL  SKIN  FRICTION 

\ 

The  similarity  laws  not  only  specify  the  velocity  profile  but  also  the  local  skin  friction. 
From  the  logarithmic  forms  of  the  inner  and  outer  similarity  laws.  Equations  (27)  and  (28) 


o = Alny  + B|  + B 


(60) 


The  objective  now  is  to  determine  a as  a function  of  local  Reynolds  number  Rg  and  shape 
parameter  H. 

Equation  (60)  is  rewritten  as 


Since  from  Equation  (34) 


and  from  Equation  (44) 


e _ 

b~  o H 


b,  = :.,-Ua 

o = AlnHR0  + B,--^A  + 2l,  - A In  I, 

1 1 is  a function  of  G from  Equations  (44)  and  (45).  An  empirical  fit  is  given  by 
21,  - A In  I,  = 1. 3462 G - 0.2605  A - 0.9392  A In  G 

Then  G is  a function  of  a and  H in  Equation  (39)  which  leads  to 


0.3462(3^889  - H)  g + Q 9392  A i„  0 = A In  R0  + B,  - 2.0938  A 


(61) 

(62) 

(63) 

(64) 

(65) 


0.93921 


. . f(H-  I)093 
" [ H 1.9392 


(66) 


This  is  the  implicit  relation  for  local  skin  friction  coefficient  a as  a function  of  Rg  and  H. 
The  factor  B,  is  a constant  for  smooth  surfaces  and  a variable  for  rough  surfaces  and/or  drag- 
reducing  polymers.  At  separation  o -*  «.  H -*  3.889. 

An  approximate  explicit  expression  may  be  obtained  which  is  accurate,  except  near 
separation.  A numerical  fit  gives 

21,  - A In  I,  * I.I67G  - 1.577  A (67) 

Then 

° ~(|.167  - O.I67h)<A  ln  H R0  + Bl  ” 3-4,0A)  (68) 


The  approximate  value  of  o may  be  used  to  obtain  a more  accurate  value  by  reiteration. 


If 


a = o + e 


where  o = approximate  value  of  a 
e = correction  of  a 

then  by  infinitesimal  analysis  of  Equation  (66) 


where 


_ 7-aq-flln  o 

a,  l 

a 


_ 0.3462(3.889  - H) 
® H 


& = 0.9392  A 


and 


7 = A In  Rfl  + B,  - 2.0938  A 


_Alnri!L^ 

l H'”9 


_ , ,0.9392' 
9392 


(69) 


(70) 


(71) 

(72) 

(73) 


Partial  derivatives  of  skin  friction  parameter  o required  in  the  subsequent  analysis  may 
be  obtained  from  Equation  (66) 


3 In  a 2.8885  | l.3462(H  - l)o 


fl.3462(H  - Do  + A H(H  - 1.9392)1 
3H  H(H  - I ) L (3.889  - H)o  + 2.71 29 AH  J 


(74) 


and 
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3 In  a = 2.8885  AH 

3 In  " <3.889  - H )o  + 2.7I29AH 


(75) 


EQUILIBRIUM  PRESSURE  GRADIENTS 
PRESSURE-GRADIENT  PARAMETER 


Similarity  is  maintained  in  the  velocity-defect  law  if  the  Clauser  pressure-gradient 
parameter3 


0 


i!  dp  _ 2 1 dU 

rw  ds  ' U ds 


(76) 


is  maintained  constant  with  respect  to  s or  d0/ds  = 0.  This  is  the  condition  for  an  equilibrium 
pressure  gradient:  B2  is  constant,  and,  for  the  most  part,  G is  constant. 

It  has  been  found  empirically3  that  for  equilibrium  pressure  gradients  G is  a function  of 
0 so  that 

G = aj  v/a2  + 0 - a3  (77) 

where  a{ , a2.  and  a-j  are  constants. 


SHAPE  PARAMETER  EQUATION 

The  local  skin  friction  coefficient  o is  given  by 


II 


o = A In  H Rfl  + B,  - — A + 2 I(  — A In  lj 


From  Equation  (39) 


H G A . „ „ B1  11  A 2,1  A , 

H-l~o~GnHR0  + G-  6 G + G ~ G > 


(64) 


(78) 


where  A is  in  general  a constant,  and  B(  is  a constant  for  smooth  surfaces. 

For  equilibrium  pressure  gradients,  G and  l|  are  constant  with  respect  to  s.  By 
differentiation  of  Equation  (78) 


AH(H  - I)2  1 d I" 

A(H-I)2+GhJ  ds 


for  equilibrium  pressure  gradients. 


(79) 


IS 


SHAPE  PARAMETER  DIFFERENTIAL  EQUATION 


The  entrainment  equation  Equation  (22)  is  converted  into  a differential  equation  for 
shape  parameter  (dH/ds)  as  follows.  In  terms  of  entrainment  shape  parameter,  H = (6  - 6*)/0 


or 


.S-E-nS-HfS 


dH  - d ln  R0 
2n=  E _ H0  . - 
ds  ds 


Now  for  two-parameter  velocity  profiles 

H = flH,  In  R#| 

Equation  (82)  when  expanded  into  partial  derivatives  becomes 


dH  _ /3H\dH  / 3H  \ d ln  Rfl 
ds  \dH  / ds  \3  In  Rq)  ds 


Then 


0 


dH_/aH\_  (m\  (r  SH  V 
ds  \3H/  \3H/  \ 3 ,n  R0/ 


d ln  R, 


0 


ds 


From  the  momentum  equation  Equation  (8) 
d 'n  fy)  1 p + "o 


0 dll 


ds 


PV 


Finally 


or 


where 


(80) 


(81) 


(82) 


(83) 


(84) 


(85) 


(86) 


„dH_  M£dU  + _M_  _L.E/  <EJL£>  (87) 

9 ds  M U ds  (H  + I)  02  \ an/  (H  + 1)  p\j2 
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is  a pressure-gradient  factor;  M represents  the  ability  of  the  shape  parameter  differential 
equation  to  respond  to  external  pressure  gradients. 

ENTRAINMENT  FACTOR 

For  equilibrium  pressure  gradients,  the  entrainment  factor  designated  as  “Ee”  is  obtained 
from  Equations  (72)  and  (78) 

or,  from  Equations  (76)  and  (85) 

where  Ee  is  a function  of  H and  Rg,  since  p depends  on  G in  Equation  (77),  and  G and  o are 
functions  of  H and  Rg  from  Equations  (39)  and  (66). 

EFFECT  OF  LOW  REYNOLDS  NUMBER 

For  smooth  flat  plates  (zero  pressure  gradient),  it  has  been  found9  that  the  velocity- 
defect  factor  B2,  instead  of  being  constant,  is  a function  of  local  skin-friction  coefficient  a 
for  low  values  of  a.  This  may  be  stated  as 

B2  = B2[0|  - AB|ol  (91) 

B2  is  the  part  of  B2  affected  by  pressure-gradient  parameter  p and  AB2  by  low  Reynolds 
number.  Since  G from  Equations  (38),  (40),  and  (41)  is  a function  of  B2,  then 

G = G[0J  - AG[o|  (92) 

G is  the  part  of  G affected  by  pressure-gradient  P and  AG  by  low  Reynolds  number. 


9 Granville,  P.S.,  “The  Drag  and  Turbulent  Boundary  Layer  of  Flat  Plates  at  Low  Reynolds  Numbers,' 
DTNSRDC  Report  4682  (Dec  1975). 


As  shown  by  Cebeci  and  Mosinskis,10  a low  Reynolds  number  effect  exists  also  for 
pressure  gradients.  For  a pressure  gradient  these  relations  may  be  generalized  to 

B2  = B2  (0|  - AB|o|  (93) 

and 

G = G[0]  - AG|o)  (94) 

Equation  (77)  also  becomes 

G[0)  = a,  v/a2  + 0 - a3  (95) 

INITIAL  CONDITIONS  AT  TRANSITION 

At  the  transition  point  from  laminar  to  turbulent  flow,  there  is  physically  no  change  in 
momentum  loss  since  such  a loss  would  require  an  impulsive  force.  Then  at  transition 

(R0>turb  = (R*>lam  (96) 

If  turbulence  stimulation  of  some  kind  were  to  be  used  to  hasten  transition,  then  an  additional 
momentum  loss,  expressed  as  AR#,  would  arise  as  a result  of  the  separation  drag  of  the 
stimulator.  Then  at  transition 


(R«'turbS(R0)lam  + AR0  (97) 

The  problem  of  an  initial  value  of  shape  parameter  H at  transition  is  more  difficult. 
During  transition  the  high  values  of  H in  laminar  flow  undergo  a rapid  decrease  to  the  lower 
values  of  H in  turbulent  flow.  For  pressure  gradients,  no  method  of  predicting  an  initial  value 
of  H seems  to  exist.  Instead,  a convenient  initial  value  of  H is  that  of  H0|R^|  for  zero 
pressure  gradient. 

For  zero  pressure  gradient,  from  Equation  (39) 


or  for  the  low  Reynolds  number  effect  of  Equation  (90) 


(98) 


Cebeci,  T.  and  GJ.  Mosinskis,  “Computation  of  Incompressible  Turbulent  Boundary  Layers  at  Low 
Reynolds  Numbers,"  American  Institute  of  Aeronautics  and  Astronautics  Journal,  Vol.  9 No  8 
pp.  1632-1634  (Aug  1971). 
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(99) 


Also  as  shown  in  Reference  3 


oQ  = aj  In  Rfl  + a2  (100) 

The  subscript  o refers  to  conditions  at  zero  pressure  gradient. 

NUMERICAL  EVALUATIONS 

Although  analytical  expressions  have  been  developed  for  the  relationships  needed  to 
determine  the  development  of  two-dimensional  turbulent  boundary  layers,  there  remain 
constants  whose  values  have  to  be  supplied  empirically. 

Numerical  comparisons  of  the  various  properties  of  the  velocity  profile  are  to  be  made 
fv  for  different  values  of  shape  parameter  H and  Reynolds  number  R#  which  represent  the  two- 

parameter  characterization.  The  value  of  Rg  = 10*  used  is  representative  of  laboratory 
. conditions,  while  the  value  of  Rg  = 106  is  representative  of  boundary  layers  found  in  many 

underwater  applications.  The  values  of  H range  from  1.4  representative  of  flat-plate  conditions 
to  3.  representing  a high  value  for  separation.  On  airfoils,  separation  may  occur  for  a value 
of  H from  1 .8  to  2.6. 

SIMILARITY  LAW  CONSTANTS 

Values  of  Bj  = 3.88  and  A = 2.606  have  been  found9  to  give  flat-plate  drag  results  close 
to  the  well-accepted  Schoenherr  line  at  high  Reynolds  numbers.  These  values  are  used  in  the 
following  numerical  evaluations. 

LOCAL  SKIN  FRICTION 

i 

Values  of  local  skin-friction  coefficient  rw/(pU2)  obtained  from  the  similarity-law 
solution  in  Equation  (66)  are  shown  in  Figure  2 as  a function  of  shape  parameter  H for  fixed 
values  of  R^  = I04  and  106.  There  is  the  expected  decrease  in  the  local  skin-friction 
coefficient  with  increasing  H,  representative  of  adverse  pressure  gradients,  and  with  increasing 
Rfl,  representative  of  a thickening  boundary  layer. 


A comparison  with  the  Ludwieg-Tillmann3  formula  is  also  shown 


Figure  2 — Wall  Shearing- Stress  Coefficient 


The  Ludwieg-Tillman  formula  represents  an  empirical  fit  to  data  around  Rg  = 104  and  is 
unreliable  for  R0  = I06. 

Also  shown  arc  test  data*1  at  = 10*  for  a variety  of  boundary  layers  which  are  listed 
in  Table  I.  There  is  good  agreement  between  the  test  data,  the  similarity  law  formula,  and 
the  Ludwieg-Tillmann  formula.  Unfortunately  there  do  not  seem  to  be  test  data  for  R0  = I06. 

TABLE  I - BOUNDARY-LAYER  DATA  FOR  R0  = I04 

(Obtained  by  interpolation  of  values  listed  in  Reference  1 1) 


Flow  No. 

Description  of  Boundary-Layer  Flow 

— xIO3 
pU2 

H 

G 

1 

H 

1100 

Ludwieg  and  Tillmann,  Mild  Adverse  Pressure  Gradient 

1.18 

1.397 

8.27 

7.61 

1200 

Ludwieg  and  Tillmann,  Strong  Adverse  Pressure  Gradient 

1.115 

1.426 

8.95 

7.76 

1500 

Tillmann,  Ledge  Flow 

1.445 

1.283 

5.80 

8.74 

2100 

Schubauer  and  Klebanoff 

1.265 

1.353 

7.34 

8.18 

2300 

Clauser,  Equilibrium  Moderate  Adverse  Pressure  Gradient 

0.57 

1.804 

18.67 

4.42 

2500 

Bradshaw,  Equilibrium  Mild  Adverse  Pressure  Gradient 

1.15 

1.426 

8.81 

6.55 

3500 

Newman,  Adverse  Pressure  Gradient 

0.755 

1.636 

14.15 

4.76 

4400 

0.785 

1.597 

13.34 

5.45 

4500 

Schubauer  and  Spangenberg,  Strong 
Adverse  Pressure  Gradient 

1.0375 

1.454 

9.69 

6.68 

4800 

1 

1.16 

1.395 

8.31 

7.05 

5000 

1 

Fraser,  Flow  in  Round  Diffuser 

0.23 

2.372 

38.14 

4.65 

5100 

; 

0.32 

2.221 

30.73 

5.13 

5300 

Stratford,  Nearly  Separating  Boundary  Layer 

0.164 

2.573 

47.74 

3.75 

Partial  derivatives  of  o,  with  respect  to  H and  R0  in  Equations  (74)  and  (75)  needed  in 
the  analysis,  are  shown  in  Figures  3 and  4.  In  comparison,  the  results  from  the  Ludwieg- 
Tillmann  formula  show  constant  values  for  (3  In  a)/3H  and  (3  In  o)/( 3 In  R0).  Experiment- 
ally derived  results  are  not  plotted  because  they  are  rot  sufficiently  precise  to  determine  such 
derivatives. 


1 1 “Proceeding -Computation  of  Turbulent  Boundary  Layers,  1968  AFOSR-IFP  Stanford  Conference,” 
Edited  by  Coles,  D.E.  and  E.A.  Hirst,  Vol.  2,  Thermosciences  Division,  Stanford  University  (1968). 
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Variation  of  d In  o/d  In  R« 


ROTTA  CLAUSER  SHAPE  PARAMETER 


The  large  numerical  variation  of  the  Rotta-Clauser  parameter  G with  H and  R0  from 
Equations  (39)  and  (66)  is  displayed  in  Figure  5.  At  separation  G -+  <». 

For  turbulent  boundary  layers  approaching  separation,  G undergoes  a large  variation. 

On  the  other  hand,  for  equilibrium  pressure  gradients.  G remains  constant  with  respect  to  the 
streamwise  direction. 

A comparison  is  also  shown  with  experimental  data  from  Table  I in  Figure  5. 

ENTRAINMENT  SHAPE  PARAMETER 

The  variation  of  entrainment  shape  parameter  H with  shape  parameter  H at  fixed  values 
ot  = I04  and  R#  - 106  is  shown  in  Figure  6 for  values  obtained  from  similarity-law 
formula;  see  Equation  (57).  It  is  to  be  noted  that  the  variation  with  is  greatest  at  the 
low  values  of  H.  For  comparison,  the  one-parameter  relation  of  Head  in  Equation  (51 ) and 
the  power-law  relation  in  Equation  (49)  arc  also  plotted  in  Figure  6. 

A comparison  with  plotted  experimental  data  listed  in  Table  1 for  Rfl  = I04  shows 
excellent  agreement  with  the  similarity-law  formula,  except  for  two  points  from  Fraser  Flows 
5000  and  5100.  Since  the  value  of  the  entrainment  shape  parameter  depends  on  the  very 
imprecise  experimental  determination  of  the  boundary-layer  thickness  (Equation  (46))  there 
is  a large  measure  of  uncertainty. 

Partial  derivatives  of  H with  respect  to  H and  Rg,  required  in  the  analysis  from  Equations 
(5H)  and  (59).  arc  shown  in  Figures  7 and  8.  For  comparison,  values  obtained  from  the  Head 
relation  of  Equation  (52)  and  the  power-law  relation  in  Equation  (50)  are  also  shown  in 
Figure  7.  Since  these  are  one-parameter  relationships,  there  is  no  variations  with  Rg. 

PRESSURE -GRADIENT  FACTOR 

The  variation  of  pressure-gradient  factor  M,  obtained  from  the  similarity  laws  of 
Equation  (88),  with  shape  parameter  H is  shown  in  Figure  9 for  fixed  values  of  R^  = 104 
and  I06.  Comparison  with  the  one-parameter  Head  and  power-law  relations  is  also  given.  At 
higher  values  of  H.  the  similarity  laws  show  somewhat  higher  values  of  M,  which  indicates  a 
more  rapid  response  to  pressure  gradient. 

There  is  also  plotted  a relation  for  M given  for  separating  boundary  layers  on  airfoils  by 
Garner,12  which  is  typical  of  older  purely  empirical  methods; 

,2Rotta.  J.C.,  “Turbulent  Boundary  Layers  in  Incompresaible  Flow,"  Chapter  in  Progress  in  Aeronautical 
Sciences,  Vol.  2,  Edited  by  A.  Ferri,  et  al„  Pergamon  Press,  New  York  (1962). 
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Figure  5 - Rotta-Clauser  Shape  Parameter 


Entrainment  Shape 


SHAPE  PARAMETER  (H) 
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Figure  9 - Variation  of  Pressure-Gradient  Factor 


M = e5(H-14> 


(102) 


The  cxponenial  form  leads  to  increasingly  large  values  of  M at  higher  values  of  H.  Such 
relationships  have  been  shown12  to  be  approximate  when  applied  to  a wide  variety  of 
pressure  gradients. 

(103) 

(104) 

(105) 


EFFECT  OF  LOW  REYNOLDS  NUMBER 

The  variation  of  AG  obtained9  empirically  may  be  fitted  by 

AG  = - O.I6o  + 3.92,  19.7  < a < 24.5 

Also 

AG  = 0.77,  o < 19.7 

AG  = 0.  a>  24.5 


EQUILIBRIUM  PRESSURE  GRADIENTS 

Since  G|0  = 0|  = 6.84  for  the  similarity  constants  used  here,  the  variation  of  G with  0 
for  equilibrium  pressure  gradient  becomes 

G = 6.1  s/  1.81  + 0 - 1.37  (106) 

which  represents  a slight  modification  of  the  relation  given  by  Nash.3 

ENTRAINMENT  FACTOR 

For  equilibrium  pressure  gradients,  the  variation  of  entrainment  factor  F,  with  shape 
parameter  H at  R#  = 104  and  106  has  been  computed  from  Equation  (90)  without  including 
the  effects  of  p and  q and  is  displayed  in  Figure  10.  Close  agreement,  is  shown  with  the 
results  for  = I04  obtained  from  the  graphs  of  Head  and  Patel.2  There  is  no  comparison 
for  = I06,  which  is  beyond  the  limit  of  the  graphs  of  Head  and  Patel. 

In  addition  the  one-parameter  values  of  E obtained  empirically  and  fitted  analytically 
are  also  shown 

Head3 

E = 0.0306  11.535  (H  - 0.7r2-715  + 0.3] -°'653  (107) 

Nicoll  and  Ramaprian3 

E = 0.035  >/H  - 1.25  ( 1 08) 
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Figure  10  - Entrainment  Factor 


CONDITION  OF  dH/ds  - 0 


A simple  solution  results  for  the  rather  stringent  requirement  that  the  shape  parameter  H 
dH 

stay  constant  or  = 0.  Then  from  Equation  (87) 


e_  du 

U ds 


(109) 


In  this  case  the  required  pressure  gradient  yy  yjy  is  a function  of  H and  Rg. 

Figure  22.1  of  Reference  12  shows  a plot  of  various  prediction  methods.  There  is  a 
tremendous  variation.  Figure  1 1 shows  a comparison  of  predictions  of  the  similarity  laws  of 
this  report  and  those  of  Head  with  experimental  data  for  equilibrium  pressure  gradients.  It  is 
evident  that  despite  the  wide  experimental  scatter,  the  similarity  laws  provide  reasonable 
agreement. 


ZERO  PRESSURE  GRADIENT  AT  LOW  REYNOLDS  NUMBERS 

For  the  low  range  of  Reynolds  numbers  800  < Rg  < 5000,  a fit  to  the  data  of 
Reference  9 gives  for  zero  pressure  gradient 

a = 3.22  In  Re  - 146  (110) 

Substituting  in  Equations  (99)  and  ( 100)  results  in  values  of  H0  plotted  in  Figure  12. 
which  compare  well  with  experimental  data.13  These  may  be  used  as  initial  values  of  H at 
transition  since  transition  usually  occurs  at  those  low  Reynolds  numbers. 


l3Landweber.  L..  "The  Frictional  Resistance  of  Flat  Plates  in  Zero  Pressure  Gradient,”  Transactions  of 
Society  of  Naval  Architects  and  Marine  Engineers,”  Vol.  61,  pp.  5-32  (1953). 
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SHAPE  PARAMETER  (H) 


SCHULTZ  GRUNOW  1 REFERENCE 


Figure  1 2 - Shape  Parameter  for  Zero  Pressure  Gradient 
at  Low  Reynolds  Numbers 


SUMMARY 


All  the  necessary  analytical  relations  have  been  derived  for  the  calculation  of  a two- 
dimensional  turbulent  boundary  layer  in  a pressure  gradient.  From  these  relations,  the 
variation  of  the  momentum  loss,  the  boundary  layer  thickness,  the  local  skin  friction,  and  the 
velocity  profile  may  be  calculated  as  a function  of  downstream  position.  The  method  is 
firmly  tounded  on  only  fundamentals,  the  equations  of  motion  for  the  boundary  layer,  the 
velocity  similarity  laws,  and  the  characteristics  of  equilibrium  pressure  gradients  which  arc 
well  substantiated  by  experimental  data.  The  two-parameter  characterization  should  be 
sufficient  for  most  engineering  analyses  of  the  two-dimensional  turbulent  boundary  layer  in 
pressure  gradients.  For  the  boundary  layer  close  to  separation,  adequate  formulations  for  the 
transverse  pressure  distribution  and  the  normal  turbulent  stresses  are  still  missing.  Basing  the 
equilibrium  entrainment  factor  on  H and  R0  only  without  a direct  coupling  to  the  pressure 
gradient  should  prove  adequate  for  most  pressure  gradients  found  in  engineering  applications. 
For  unusually  large  pressure  gradients,  however,  adjustments  may  be  needed  to  the  equilibrium 
entrainment  factor. 

For  convenience,  the  necessary  analytical  relations  are  listed  in  the  appendix  for  the 
solution  of  the  two  simultaneous  differential  equations:  the  momentum  and  the  shape 
parameter.  Also  included  are  the  relations  specifying  the  boundary-layer  velocity  profile  at 
each  streamwisc  position. 


APPENDIX  A 

CALCULATION  PROCEDURE 

All  the  necessary  analytical  relations  for  calculating  the  development  fo  the  two- 
dimensional  turbulent  boundary  and  the  associated  velocity  profiles  are  now  listed. 

Solve  as  simultaneous  equations  for  0[s)  and  H[s],  given  U(s|  and  ^(s)  <P  and  o 
are  assumed  equal  to  zero). 


d6 

ds 


(H  + 2)£  j!  + J- 

U ds  02 


(4) 


where 


(87) 


(88) 

(90) 

(95) 


G = G + AG  (94) 

AG  = 0.16  a + 3.92,  19.7  < a < 24.5  (102) 

AG  = 0.77,  a < 19.7  (103) 

AG  = 0.  a > 24.5  (104) 

g = ° _ n)  (39) 


~ ( H2  \ / 0.4739  A 5 A2  \ „ 

H V 4857  + ~ G~  + tfTs)- » 

aFI  H(H  - 2)  - H / H2  \ / 0.4739  A 13.75  A2\/  1 d In  o\ 

3H  H(H-I)  ‘^H-iy^  G + c2.75  y\H(H-l)  + 3H  J 


(58) 


3H 

3 In  Rfl 


(59) 


/ H2  \ / 0.4739  A 13.75  aA  3 In  q 

G + C2.7S  jdlnRfl 

0.3462(^.889  ^_Hj  o + o qv)2  A |n  0 = A In  R0  + B,  - 2.0938  A - A In 


) In  q _ 18885  p.3462(H  - l)q  + A H(H  - 1.9392)1 
3H  H<  H - 1 >L  (3.889  - H)q  + 2.7129  AH  J 


3 In  q 


2.8885  AH 


3 In  Re  (3.889  - H)q  + 2.7129  AH 

A = 2.606.  B,  = 3.88 
a(  = 6. 1.  a2  = 1.81,  a^  = 1.37 


(66) 

(74) 

(75) 


INITIAL  CONDITIONS  AT  TRANSITION 


(U0\ 


turb 


(R0*,urb  = (Re>.»m  + ^Ri 


lam 


I G0-AG|q0] 

H0"  ' 


G0  = 6.84 


AG  = - 0.l6qo  + 3.92 


oG  = 3.22  In  Rfl  - 1.46 


s e 

VELOCITY  PROFILE 

At  each  station  s,  H.  9.  a and  G have  been  calculated.  Then  from  Equation  (56) 

G 


(97) 

(99) 

(102) 

(109) 


»1  " 


from  Equation  (44) 


G q2.75 


B2  = 2'l 
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From  definitions 


From  Equation  (60) 


Then 


For  y*  > y0* 


— = A In  y*  + B,  + B,  ^ + Aq, 


For  y*  < y * ( Reference  6) 

laminar  sublayer,  = y*  o < y*  < y^ 

buffer  sublayer  — = A In  (y*  - J)  + Bj , < y*  < 

UT 

J = Bj  + A In  A - A 
y£  = B,  + A In  A 

More  simply  let 

— = A In  (y*  - J)  + B,  + B2  -5-  + Aq,  yL*  < y*  < tj 
uT  L 

for  the  boundary  layer  outside  the  laminar  sublayer. 
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